3.3.28 \(\int \frac {(f x)^m (d+e x^2)}{(a+b x^2+c x^4)^{3/2}} \, dx\) [228]

3.3.28.1 Optimal result
3.3.28.2 Mathematica [A] (verified)
3.3.28.3 Rubi [A] (verified)
3.3.28.4 Maple [F]
3.3.28.5 Fricas [F]
3.3.28.6 Sympy [F]
3.3.28.7 Maxima [F]
3.3.28.8 Giac [F]
3.3.28.9 Mupad [F(-1)]

3.3.28.1 Optimal result

Integrand size = 29, antiderivative size = 323 \[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\frac {d (f x)^{1+m} \sqrt {1+\frac {2 c x^2}{b-\sqrt {b^2-4 a c}}} \sqrt {1+\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}} \operatorname {AppellF1}\left (\frac {1+m}{2},\frac {3}{2},\frac {3}{2},\frac {3+m}{2},-\frac {2 c x^2}{b-\sqrt {b^2-4 a c}},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}\right )}{a f (1+m) \sqrt {a+b x^2+c x^4}}+\frac {e (f x)^{3+m} \sqrt {1+\frac {2 c x^2}{b-\sqrt {b^2-4 a c}}} \sqrt {1+\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}} \operatorname {AppellF1}\left (\frac {3+m}{2},\frac {3}{2},\frac {3}{2},\frac {5+m}{2},-\frac {2 c x^2}{b-\sqrt {b^2-4 a c}},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}\right )}{a f^3 (3+m) \sqrt {a+b x^2+c x^4}} \]

output
d*(f*x)^(1+m)*AppellF1(1/2+1/2*m,3/2,3/2,3/2+1/2*m,-2*c*x^2/(b-(-4*a*c+b^2 
)^(1/2)),-2*c*x^2/(b+(-4*a*c+b^2)^(1/2)))*(1+2*c*x^2/(b-(-4*a*c+b^2)^(1/2) 
))^(1/2)*(1+2*c*x^2/(b+(-4*a*c+b^2)^(1/2)))^(1/2)/a/f/(1+m)/(c*x^4+b*x^2+a 
)^(1/2)+e*(f*x)^(3+m)*AppellF1(3/2+1/2*m,3/2,3/2,5/2+1/2*m,-2*c*x^2/(b-(-4 
*a*c+b^2)^(1/2)),-2*c*x^2/(b+(-4*a*c+b^2)^(1/2)))*(1+2*c*x^2/(b-(-4*a*c+b^ 
2)^(1/2)))^(1/2)*(1+2*c*x^2/(b+(-4*a*c+b^2)^(1/2)))^(1/2)/a/f^3/(3+m)/(c*x 
^4+b*x^2+a)^(1/2)
 
3.3.28.2 Mathematica [A] (verified)

Time = 11.34 (sec) , antiderivative size = 307, normalized size of antiderivative = 0.95 \[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\frac {x (f x)^m \left (-b+\sqrt {b^2-4 a c}-2 c x^2\right ) \sqrt {\frac {b-\sqrt {b^2-4 a c}+2 c x^2}{b-\sqrt {b^2-4 a c}}} \left (\frac {b+\sqrt {b^2-4 a c}+2 c x^2}{b+\sqrt {b^2-4 a c}}\right )^{3/2} \left (d (3+m) \operatorname {AppellF1}\left (\frac {1+m}{2},\frac {3}{2},\frac {3}{2},\frac {3+m}{2},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}},\frac {2 c x^2}{-b+\sqrt {b^2-4 a c}}\right )+e (1+m) x^2 \operatorname {AppellF1}\left (\frac {3+m}{2},\frac {3}{2},\frac {3}{2},\frac {5+m}{2},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}},\frac {2 c x^2}{-b+\sqrt {b^2-4 a c}}\right )\right )}{\left (-b+\sqrt {b^2-4 a c}\right ) (1+m) (3+m) \left (a+b x^2+c x^4\right )^{3/2}} \]

input
Integrate[((f*x)^m*(d + e*x^2))/(a + b*x^2 + c*x^4)^(3/2),x]
 
output
(x*(f*x)^m*(-b + Sqrt[b^2 - 4*a*c] - 2*c*x^2)*Sqrt[(b - Sqrt[b^2 - 4*a*c] 
+ 2*c*x^2)/(b - Sqrt[b^2 - 4*a*c])]*((b + Sqrt[b^2 - 4*a*c] + 2*c*x^2)/(b 
+ Sqrt[b^2 - 4*a*c]))^(3/2)*(d*(3 + m)*AppellF1[(1 + m)/2, 3/2, 3/2, (3 + 
m)/2, (-2*c*x^2)/(b + Sqrt[b^2 - 4*a*c]), (2*c*x^2)/(-b + Sqrt[b^2 - 4*a*c 
])] + e*(1 + m)*x^2*AppellF1[(3 + m)/2, 3/2, 3/2, (5 + m)/2, (-2*c*x^2)/(b 
 + Sqrt[b^2 - 4*a*c]), (2*c*x^2)/(-b + Sqrt[b^2 - 4*a*c])]))/((-b + Sqrt[b 
^2 - 4*a*c])*(1 + m)*(3 + m)*(a + b*x^2 + c*x^4)^(3/2))
 
3.3.28.3 Rubi [A] (verified)

Time = 0.53 (sec) , antiderivative size = 323, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.069, Rules used = {1674, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\left (d+e x^2\right ) (f x)^m}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx\)

\(\Big \downarrow \) 1674

\(\displaystyle \int \left (\frac {d (f x)^m}{\left (a+b x^2+c x^4\right )^{3/2}}+\frac {e (f x)^{m+2}}{f^2 \left (a+b x^2+c x^4\right )^{3/2}}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {d (f x)^{m+1} \sqrt {\frac {2 c x^2}{b-\sqrt {b^2-4 a c}}+1} \sqrt {\frac {2 c x^2}{\sqrt {b^2-4 a c}+b}+1} \operatorname {AppellF1}\left (\frac {m+1}{2},\frac {3}{2},\frac {3}{2},\frac {m+3}{2},-\frac {2 c x^2}{b-\sqrt {b^2-4 a c}},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}\right )}{a f (m+1) \sqrt {a+b x^2+c x^4}}+\frac {e (f x)^{m+3} \sqrt {\frac {2 c x^2}{b-\sqrt {b^2-4 a c}}+1} \sqrt {\frac {2 c x^2}{\sqrt {b^2-4 a c}+b}+1} \operatorname {AppellF1}\left (\frac {m+3}{2},\frac {3}{2},\frac {3}{2},\frac {m+5}{2},-\frac {2 c x^2}{b-\sqrt {b^2-4 a c}},-\frac {2 c x^2}{b+\sqrt {b^2-4 a c}}\right )}{a f^3 (m+3) \sqrt {a+b x^2+c x^4}}\)

input
Int[((f*x)^m*(d + e*x^2))/(a + b*x^2 + c*x^4)^(3/2),x]
 
output
(d*(f*x)^(1 + m)*Sqrt[1 + (2*c*x^2)/(b - Sqrt[b^2 - 4*a*c])]*Sqrt[1 + (2*c 
*x^2)/(b + Sqrt[b^2 - 4*a*c])]*AppellF1[(1 + m)/2, 3/2, 3/2, (3 + m)/2, (- 
2*c*x^2)/(b - Sqrt[b^2 - 4*a*c]), (-2*c*x^2)/(b + Sqrt[b^2 - 4*a*c])])/(a* 
f*(1 + m)*Sqrt[a + b*x^2 + c*x^4]) + (e*(f*x)^(3 + m)*Sqrt[1 + (2*c*x^2)/( 
b - Sqrt[b^2 - 4*a*c])]*Sqrt[1 + (2*c*x^2)/(b + Sqrt[b^2 - 4*a*c])]*Appell 
F1[(3 + m)/2, 3/2, 3/2, (5 + m)/2, (-2*c*x^2)/(b - Sqrt[b^2 - 4*a*c]), (-2 
*c*x^2)/(b + Sqrt[b^2 - 4*a*c])])/(a*f^3*(3 + m)*Sqrt[a + b*x^2 + c*x^4])
 

3.3.28.3.1 Defintions of rubi rules used

rule 1674
Int[((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + ( 
c_.)*(x_)^4)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(f*x)^m*(d + e*x^2)^q* 
(a + b*x^2 + c*x^4)^p, x], x] /; FreeQ[{a, b, c, d, e, f, m, p, q}, x] && N 
eQ[b^2 - 4*a*c, 0] && (IGtQ[p, 0] || IGtQ[q, 0] || IntegersQ[m, q])
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.3.28.4 Maple [F]

\[\int \frac {\left (f x \right )^{m} \left (e \,x^{2}+d \right )}{\left (c \,x^{4}+b \,x^{2}+a \right )^{\frac {3}{2}}}d x\]

input
int((f*x)^m*(e*x^2+d)/(c*x^4+b*x^2+a)^(3/2),x)
 
output
int((f*x)^m*(e*x^2+d)/(c*x^4+b*x^2+a)^(3/2),x)
 
3.3.28.5 Fricas [F]

\[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\int { \frac {{\left (e x^{2} + d\right )} \left (f x\right )^{m}}{{\left (c x^{4} + b x^{2} + a\right )}^{\frac {3}{2}}} \,d x } \]

input
integrate((f*x)^m*(e*x^2+d)/(c*x^4+b*x^2+a)^(3/2),x, algorithm="fricas")
 
output
integral(sqrt(c*x^4 + b*x^2 + a)*(e*x^2 + d)*(f*x)^m/(c^2*x^8 + 2*b*c*x^6 
+ (b^2 + 2*a*c)*x^4 + 2*a*b*x^2 + a^2), x)
 
3.3.28.6 Sympy [F]

\[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\int \frac {\left (f x\right )^{m} \left (d + e x^{2}\right )}{\left (a + b x^{2} + c x^{4}\right )^{\frac {3}{2}}}\, dx \]

input
integrate((f*x)**m*(e*x**2+d)/(c*x**4+b*x**2+a)**(3/2),x)
 
output
Integral((f*x)**m*(d + e*x**2)/(a + b*x**2 + c*x**4)**(3/2), x)
 
3.3.28.7 Maxima [F]

\[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\int { \frac {{\left (e x^{2} + d\right )} \left (f x\right )^{m}}{{\left (c x^{4} + b x^{2} + a\right )}^{\frac {3}{2}}} \,d x } \]

input
integrate((f*x)^m*(e*x^2+d)/(c*x^4+b*x^2+a)^(3/2),x, algorithm="maxima")
 
output
integrate((e*x^2 + d)*(f*x)^m/(c*x^4 + b*x^2 + a)^(3/2), x)
 
3.3.28.8 Giac [F]

\[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\int { \frac {{\left (e x^{2} + d\right )} \left (f x\right )^{m}}{{\left (c x^{4} + b x^{2} + a\right )}^{\frac {3}{2}}} \,d x } \]

input
integrate((f*x)^m*(e*x^2+d)/(c*x^4+b*x^2+a)^(3/2),x, algorithm="giac")
 
output
integrate((e*x^2 + d)*(f*x)^m/(c*x^4 + b*x^2 + a)^(3/2), x)
 
3.3.28.9 Mupad [F(-1)]

Timed out. \[ \int \frac {(f x)^m \left (d+e x^2\right )}{\left (a+b x^2+c x^4\right )^{3/2}} \, dx=\int \frac {{\left (f\,x\right )}^m\,\left (e\,x^2+d\right )}{{\left (c\,x^4+b\,x^2+a\right )}^{3/2}} \,d x \]

input
int(((f*x)^m*(d + e*x^2))/(a + b*x^2 + c*x^4)^(3/2),x)
 
output
int(((f*x)^m*(d + e*x^2))/(a + b*x^2 + c*x^4)^(3/2), x)